A transactional theory of fluctuations in company size 
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Detailed empirical studies of publicly traded business firms have established that the standard 
deviation of annual sales growth rates decreases with increasing firm sales as a power law, and that 
the sales growth distribution is non-Gaussian with slowly decaying tails. To explain these empirical 
facts, a theory is developed that incorporates both the fluctuations of a single firm's sales and the 
statistical differences among many firms. The theory reproduces both the scaling in the standard 
deviation and the non-Gaussian distribution of growth rates. Earlier models reproduce the same 
empirical features by splitting firms into somewhat ambiguous subunits; by decomposing total sales 
into individual transactions, this ambiguity is removed. The theory yields verifiable predictions and 
accommodates any form of business organization within a firm. Furthermore, because transactions 
are fundamental to economic activity at all scales, the theory can be extended to all levels of the 
economy, from individual products to multinational corporations. 

PACS numbers: 89.65.Gh,89.75.Da,87.23.-n 



I. INTRODUCTION 

In 1931, Gibrat introduced a formal model of market 
structure based on his Law of Proportionate Effect, which 
describes the time evolution of business firm size [1, 2]. 
There is no unique measure of the size of a company; 
a firm's size is frequently defined to be the total an- 
nual sales. Gibrat postulated that all firms, regardless 
of size, have an equal probability to grow a fixed frac- 
tion in a given period, i.e., the firm size undergoes a 
multiplicative random walk. Over appropriately chosen 
periods, the distribution of the annual fractional size 
change is approximately Gaussian, independent of the 
firm size [2] . Gibrat hypothesized that the distribution of 
firm sizes follows a nonstationary lognormal distribution. 
Subsequent statistical studies [3, 4] have found that the 
firm-size distribution is well-approximated by a lognor- 
mal distribution as predicted. More recent studies [5-10] 
have also demonstrated three empirical facts incompati- 
ble with Gibrat 's hypothesis: (i) Firm growth rates, de- 
fined to be the annual change in the logarithm of sales, 
follow a non-Gaussian distribution with slowly decaying 
tails, (ii) The standard deviation of growth rates scales 
as a power law with firm sales, e.g., a g ~ >S~' 9 , where 
/3 ~ 0.2 [5], (iii) The firm size distribution is approxi- 
mately stationary over tens of years. 

A number of models have been proposed to explain 
these incompatibilities between Gibrat 's theory and the 
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empirical facts. These models decompose each firm into 
fluctuating subunits [11-19], which are usually taken to 
be management divisions or product lines. While the 
subunits obey simple dynamics, the aggregate entities 
approximate the empirically observed statistics. Studies 
ranging from the gross domestic product (GDP) of coun- 
tries to the sales of individual products have found com- 
parable scaling laws and non-Gaussian distributions [16- 
21]. While the statistics of GDP growth may be ex- 
plained in terms of the contributions of individual firms, 
the firm sales do not obey simple dynamics. Likewise, 
while firm sales growth may be explained in terms of 
the contributions from products, products exhibit com- 
plex behavior similar to the entire firm [16]. Empirically, 
no economic entity obeying simple dynamics has been 
identified. In order to reconcile the empirical results 
of Refs. [16-21] with these theoretical models, either a 
subunit obeying a postulated simple dynamics must be 
identified, or we must account for the complex nature of 
the subunits themselves, as illustrated by the hierarchal 
schematic in Fig. 1. 

Here we develop a theory to explain the observed firm 
growth statistics that does not invoke an inherently hi- 
erarchical model of the economy [5, 11, 35]. We do not 
make any assumptions about the internal structure of 
firms. We postulate that the total annual sales of a par- 
ticular firm can be broken down into a finite sum of indi- 
vidual sales transactions. In our analysis, each firm's to- 
tal annual sales is an independent random variable that 
is characterized by three parameters (see below). Fur- 
thermore, these parameters vary randomly from firm to 
firm. We refer to the set of firms as the population. We 



2 



(^Country 2 




FIG. 1: The GDP of a country is decomposed into the contri- 
butions of individual firms. Firms' contributions are further 
split by individual products, ad infinitum. Because complex 
behavior is observed at all scales, models of economies that ex- 
plain complex behavior in terms of fluctuating subunits must 
account for the composite nature of the subunits themselves. 



define heterogeneity as a measure of the variability of the 
firms' parameters within a population. 

In Section II, we develop a model of a single firm's 
sales. In Section III, we quantify heterogeneity in a pop- 
ulation of firms and derive a scaling relationship between 
firm size and the standard deviation of growth rates. Sec- 
tion IV studies the distribution of annual growth rates 
for a single firm and approximates the growth distribu- 
tion generated by a heterogeneous population. Section V 
uses the scaling analysis presented in Sec. Ill to illustrate 
the relationship between the number of products sold by 
firms and firm sales. 



II. MODEL FOR THE STATISTICS OF A 
SINGLE FIRM 

We define the size of firm i to be its total annual sales 
Si which is the sum of Ni transactions, 



Si = Tii + Ti 



i,Ni, 



(1) 



where Si > and T^k > are measured in units of 
currency. We assume both Ni and T^ are independent 
random variables with finite moments. The total sales Si 
is then a random variable. Firm i may increase its total 
sales by increasing its sales volume Ni, or by selling more 
expensive products (increasing T). The uncertainties in 
price and number depend on the nature of the products 
that a firm sells and the market. We consider Poisson- 
distributed demand [22-25], where Ni is described by a 
Poisson distribution, and we define 



Qi = (Ni) = var[7V 4 ] . 



(2) 



With this choice, the size of a firm Si is a random vari- 
able with a compound Poisson distribution, the details 
of which depend on the statistics of the individual trans- 
actions T. 



For a single firm, the time-averaged mean and variance 
of sales are denoted by 



(Si) — Qi Xi 
var[5i] = Qi Xfa. 



(3) 



where Xi = (T) is the mean transaction size and 
Cj = 1 + var[Tj] /Xf is a measure of the statistical dis- 
persion in transaction sizes. Because of its importance 
in the subsequent analysis, we introduce the logarithm 
of the firm size, 



Si = log Si. 



(4) 



The log size s, is a random variable that is a function 
of the random variable Si. To simplify the notation, we 
also introduce the logarithms of the expected number of 
transactions 



qi = \ogQ t , 
and the mean transaction size 

Xi = \ogXi. 



(5) 



(6) 



We estimate the first moment of Sj by computing the 
expected value of the Taylor series expansion of Eq. (4) 
about the mean (Si), 



(si) = log (QiXi)-^- + 



(7) 



In the limit that the number of transactions is large, 
Qi 3> 1, we retain only the leading order term, 



(Si) fts log (Qi Xi) = qi+Xi. 
The variance of is estimated similarly, 



var[si] 



Q t (Xi^Cj 

QUXi) 2 



Q t 



(8) 



(9) 



Importantly, the variance in the log size of a firm is in- 
versely proportional to the number of transactions. We 
combine Eqs. (4), (8), and (9), to compactly express the 
(approximate) log size of a firm in terms of its mean and 
standard deviation, 



Si = (Si) + (var[sj] ) 1/2 r) 

1/2 



= q t +Xi + 



Qi 



(10) 



where fj is a random variable that has zero mean and 
unit variance. Equation (10) expresses the log sales of a 
particular firm described by the three parameters qi, x i} 
and Cj. 

To compute the statistics of sales growth, recall our 
assumption that the annual sales of firm i is statistically 
independent of its sales in the previous year. Let s, and 
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denote the log size of a firm in two consecutive years, 
then the annual growth is 



9i = s i~ s i 



1/2 



{ff - if) . 



(11) 



The growth gi is a random variable with vanishing odd- 
integer moments and variance given by 



2c- 

var [ft] = 7T = 2c i exp(-%). 



(12) 



Thus, we have defined a firm's total sales to be a sum 
of a random number of transactions of random size. Each 
firm is uniquely described by its average sales volume and 
its transaction sizes, characterized by qi, Xi, and Cj. The 
statistical independence of sales from year to year implies 
that the variance in growth rates, Eq. (12), is a simple 
function of average sales volume and the transaction size 
dispersion. 



III. QUANTIFYING HETEROGENEITY IN 
THE POPULATION OF FIRMS 



Prior studies [5-10] examined the statistics of sales 
growth within a population of firms and found that the 
standard deviation of growth rates scales as a power law 
with firm sales, a g ~ S^ 13 , and that the distribution of 
growth rates is non-Gaussian with slowly decaying tails. 
Analogously, we examine the statistics of sales growth for 
a population of firms. We will find that the heterogene- 
ity of firms is fundamental to reproducing the power-law 
scaling in the standard deviation of growth rates. 

We assume that an individual firm's parameters, q%,Xi, 
and Ci, evolve slowly and may be treated as fixed quan- 
tities. To introduce heterogeneity within the population, 
these parameters are sampled from a time-independent 
distribution [28]. In the subsequent discussion, we drop 
the subscript i if a parameter is a random variable that 
is representative of all firms (not one firm in particu- 
lar.) Our analysis assumes that the parameter c has finite 
mean (c), and that it is statistically independent of the 
other firm parameters. The means and variances of the 
parameters q and x across the population are denoted by 



(q) = TYlq 

(x) = m x 



var[ 9 ] = a\ 
var[x] = a 2 x . 



(13) 
(14) 



We take both q and x to be normally distributed within 
the population, with correlation coefficient p. The joint 
probability density function (PDF) of the pair is given 
by 



p(q, x) 



exp [— 2^(1 



4ttW(1-P 2 ) 



(15) 




x, log transaction size 

FIG. 2: Solid ellipses are the level curves of the bivariate dis- 
tribution given in Eq. (15) with p = 0. Dashed lines show pos- 
sible configurations of transaction volume and size for firms 
of various sizes. The arrow indicates that the mean log trans- 
action quantity increases linearly with log size, in agreement 
with Eq. (19). 



where 



_ (q - m q ) 2 (x - m x f 

Z — o T o 



2p(q 



m,) (x 



(16) 



<?a Co- 



The correlation coefficient depends on the industry: a 
negative value of p implies that firms that sell smaller 
quantities tend to sell more expensive items. The pa- 
rameters of Eq. (15) also reflect the fact that firms must 
minimize inefficiency and risk. Firms with too many 
transactions have excess overhead, and firms with low 
volume risk bankruptcy from large sales fluctuations. 

We examine the statistics of firms of a fixed mean size, 
s av g = q + x = const. As illustrated in Fig. 2, the rela- 
tionship between total firm sales and transactions imply 
that statistics that depend on transaction volume qi also 
depend on the sizes of the firms under consideration. For 
example, to extend the variance in growth rates given 
in Eq. (9) to include all firms of a fixed size, we must 
determine the distribution of q conditioned on s avg , 



p(q\s. 



avgy 



P(q, s avg) 

J ' dqpiq.s^Y 



(17) 



where p(q, s a vg) is a Gaussian PDF that is obtained by 
substituting the constraint, x = s avg — q, into Eq. (15), 



p(<7,Sav g ) =p(q,x = s 



avg 



(18) 



The quotient of two Gaussian distributions is itself 
Gaussian. Here, the conditional distribution in Eq. (17) 
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is completely characterized by its conditional mean 

avg 



+ 



(m q - m x )pa q a x + m q a 2 - m x a 2 



and conditional variance 
var[g|s avg ] 



(1 - p 2 )a 2 a x 



where 



a + pa q a x 
0= , u = a q 



2pa q a x . 



(19) 



(20) 



(21) 



We use Eqs. (12), (19), and (20) to find the variance 
in annual growth rates averaged over all firms of a given 
log size, 



var[g | 



2(c) exp 



I s avg ) + 



v&r{q | s 



• (22) 



This result recovers the empirically observed power-law 
relationship between firm sales S avg = exp(s avg ) and vari- 
ability in growth rates, 



var^lSavg] ~ S 2(} 



(23) 



Among all firms of a fixed mean size, we find that the 
mean number of transactions per year scales as a power 
law, 



(Q\s, 



avg/ 



a 2/3 
^avgi 



(24) 



where Q = e q . A similar calculation for the distribution 
of the mean transaction size conditioned on firm size also 
gives a power law, 



(X | Savg) - S^ 20 , 



(25) 



where X = e x . 

Finally, we note that empirical data arc frequently 
sparse and are binned into large size classes. Appendix A 
reviews how binning affects the conditional statistics pre- 
sented in this section. 



IV. DISTRIBUTION OF GROWTH RATES 

We next address the empirically observed non- 
Gaussian distribution of sales growth rates [5-10]. In 
general, a firm's sales growth distribution p gi depends on 
a set of firm-specific parameters, which we denote with 
the vector a,. For a set of heterogeneous firms T, let 
p a (S) denote the joint PDF of the parameters. Schemat- 
ically, the growth rate distribution computed for all firms 
in the set T is given by 



^■{g\^) = j da,p gl {g\a)p a {a\T) 

P 00 

= / d* 2 p(g\a 2 )p(a 2 \f), 
Jo 



(26) 



where the second equality is obtained by integrating over 
all degrees of freedom except the variance in growth 
rates, a 2 . For a sufficiently broad variance distribu- 
tion p(a 2 \!F), the distribution of growth rates II (5 is 
non-Gaussian with slowly decaying tails, a typical fea- 
ture of growth models that incorporate heterogeneity of 
firms [12, 14, 15, 17-19, 30]. 

In the present treatment, we note that within a pop- 
ulation of firms of fixed average annual sales, there are 
firms with large numbers of small transactions and vice 
versa. Because the variance in growth rates is inversely 
proportional to the number of transactions, the hetero- 
geneity in transaction volume gives rise to a distribution 
in the variance of firm growth rates. To study the dis- 
tribution of growth rates beyond the second moment, we 
must estimate the distribution of a firm's total sales Si. 

In actuarial science and operations research litera- 
ture [25, 26], the distribution of Si has been approxi- 
mated by the truncated normal, lognormal, and gamma 
distributions. For analytical convenience, we approxi- 
mate the compound Poisson distribution of firm sales Si 
with a gamma distribution. As required by our theory, 
the gamma distribution admits only nonnegative sales, 
Si > 0, and it is asymptotically Gaussian when the num- 
ber of transactions is large, in agreement with the central 
limit theorem. The PDF of firm size & is assumed to be 



PSi 



sr 



bV r(o0 



exp 



Si 
b, 



(27) 



where the parameters and bi are unique to firm i. 
These parameters are determined by equating the mean 
and variance to their values given in Eq. (3), 



(Si) = QiXi = b t ai, 
var[Si] = QiX 2 a = b 2 ^. 



(28) 



We solve Eq. (28) to obtain an expression for at in terms 
of the firm parameters, 



Qi _ exp(%) 



The PDF of a firm's log size s» follows from Eq. (27), 
1 



(29) 



Psi — , ai 



exp[ajSi - exp(sj)/6j]. (30) 



Equation (11) defines the growth of a single firm to be the 
difference of two independent random variables. There- 
fore, the PDF of single-firm growth rates is given by 



dg' Psi(g') Psi(g' + g t ), 



r(2a. t + 1) 
a2 a + 1 T 2 (a i ) 



(1 + cosh a% 



(31) 



£(^ + |) 2a 

2V^r(o 4 ) 



■(f)- 



5 




24681012 -2 -1 1 2 



s avg (log sales) g (sales growth) 

FIG. 3: Simulation of 3000 firms for 60 periods with m q = 4.0, cr, = 2.0, m x — 3.0, = 3.0, p = 0, and constant c = 1.01. (a) 
Power-law scaling of the standard deviation of growth rates conditioned on size, and fit to Eq. (23). (b) PDFs of growth rates 
for firms of three sizes and the corresponding approximate PDFs from Eq. (35). 



In the case considered here, the distribution of single firm 
growth rates depends on a single scalar parameter, a%. 
The variance in a firm's growth rate is given by 

2 1 

var^] =(gf) = 2*i(o<) + — + •••, (32) 

a>i a.i 

where Vti is the first derivative of the digamma func- 
tion [27], 

* x {x) = ^{x) = ^\ogT{x), (33) 
ax ax z 

In the large-ai limit, we recover the expression for the 
variance in Eq. (12). This asymptotic agreement is ex- 
pected because in Sec. II we employed a low-order ex- 
pansion of the logarithm. Furthermore, because the hy- 
perbolic cosine is asymptotically exponential, the tails of 
the distribution in Eq. (31) decay exponentially, 

The distribution of growth rates for a population must 
reflect the variability of the parameter a among the firms. 
Formally, the PDF of growth rates for a population is 
given by a weighted mixture of the PDFs of single-firm 
growth rates given in Eq. (31). To express this distri- 
bution analytically, we assume negligible variance in the 
transaction size dispersion, i.e., a = c is constant. In this 
case, the distribution of the parameter a is lognormal, 



1 



Pa = 



exp 



- — (logo - M) 2 
2V K ' 



(34) 



where M = {q\T) — logc and V = var^l^ 7 ] are deter- 
mined by the statistics of J 7 , the set of firms under con- 
sideration. For example, to determine the growth rate 
distribution for the entire population, we use the mo- 
ments defined in Eq. (13). Alternatively, to compute the 
distribution of growth rates among firms of fixed size, we 
use the conditional moments given in Eqs. (19) and (20). 



In the limit of many transactions, a» 1, we approximate 
the PDF of growth rates for a set of heterogeneous firms, 

exp \UVM + V 2 - Aw 2 - 8w) 1 8Vl 



where 



w 



2V 



log (cosh | D ■ 1 



(36) 



and W is the Lambert W function [29], defined by 
W(x)e w ^ = x. We give the details of the derivation 
and an asymptotic analysis of Eq. (35) in Appendix B. 

Figure 3 compares Eq. (23) and Eq. (35) to the statis- 
tics generated from a simulation of individual firms with 
annual sales sampled according to Eq. (27). We see that 
our analysis produces a tent-shaped distribution, consis- 
tent with empirical facts [5, 6, 9]. 



V. THE SCALING OF PRODUCTS 

Power-law scaling relationships are observed in aca- 
demic, ecological, technological, and economic systems 
that relate a measure of size to the number of con- 
stituents [31-37]. For example, the number of dif- 
ferent products sold by firms (the constituents) grows 
as a power law in firm sales (a measure of size), 
AT~ S QA2 [16]. 

To elucidate the relationship between products and 
sales, we adopt the notation of Section III. Let Qi de- 
note the number of products within firm i and let Xi de- 
note the firm's total annual sales per product. Because 
economic conditions are different in different locales and 
industries, we assume that among all firms, the number 
of products Q and the mean sales per product X are 
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both lognormally distributed random variables. For con- 
venience, we drop the subscript on the random variables 
and take the logarithm, 

S = QX -> s = q + x. (37) 

The expression for the joint density of q and x is identical 
to Eq. (15), except that in this context, the parameters of 
the distribution reflect aspects of the product portfolios 
of firms. 

From Eqs. (24) and (25), we find that the mean number 
of products and the mean product size scale as power laws 
with average firm size, 

(Q|S avg > ~ (38) 
(X\S avg ) ~ 5 a 1 " 2 ' 3 , (39) 

where /3 is given in Eq. (21). With a suitable choice of 
the parameters cr 2 , cr 2 and p, one can reproduce the ob- 
served power-law scaling relationship between the num- 
ber of products and firm sales. Moreover, for p = and 
with cr 2 — W, and cr 2 = V, the analysis here is identical 
to that of Refs. [12, 16]. 

VI. DISCUSSION 

The present model of firm size fluctuations produces a 
non-Gaussian distribution of growth rates and generates 
a stationary lognormal distribution of firm sizes, features 
consistent with empirical studies. Our approach bears 
similarity to the model proposed in Ref. [12] in which a 
firm is split into subunits that each obey simple dynam- 
ics. In the present work, the number of transactions is 
analogous to the number of subunits, and the transaction 
size parallels the sizes of the subunits within a firm. The 
model of Ref. [12] postulates a complex internal structure 
within a firm. In contrast, our approach decomposes a 
firm's sales into the individual transactions that occur in 
a year. 

Because transactions are well-defined measurable 
quantities, the assumptions and predictions of the theory 
are verifiable. In practice, correlations between transac- 
tions invalidate our assumption that the number of trans- 
actions Ni follows a Poisson distribution. Indeed, trans- 
actions within a firm can be regrouped to remove cor- 
relations. Because these uncorrelatcd groups satisfy the 
assumptions of our model, we conclude that a Poisson- 
distributed number of transactions is not strictly neces- 
sary to obtain the results presented here. To retain our 
transactional framework, one can extend the model of a 
single firm's statistics described in Section II to account 
for correlations between transactions at the expense of 
additional parameters [24] . Independent of the statistical 
details, most measures of microeconomic and macroeco- 
nomic consumption tally individual transactions and are 
therefore amenable to the transactional approach pre- 
sented here. Figuratively speaking, transactions are the 
'atoms' of economic activity: they are discrete and are 



the basis of economies of all scales, from the individual 
to the national. 

Using our approach, we studied the relationship be- 
tween the number of different products sold by firms and 
total sales receipts. We found that differences between 
firms could account for the reported scaling relationship 
between the two quantities, (Q |5 avg ) ~ S&vg- O ur anal- 
ysis applies to many systems where two distinct mea- 
sures of size are related by a proportionality constant that 
varies because of economic, environmental, and other ac- 
cidental circumstances. When the relationship between 
the measures is weak, the population is heterogeneous 
and non-trivial power laws emerge [35]. For example, 
both transaction volume and total sales quantify the size 
of a firm, but no strict universal proportionality between 
these quantities exists. In this case, Eq. (24) predicts a 
power-law relationship between the quantities. In ecol- 
ogy, both island surface area and number of species quan- 
tify island size: islands with greater area generally har- 
bor more species. However, the mix of species and the 
number distinct niches varies from island to island. In 
agreement with our hypothesis, empirical studies show 
that the number of species scales as a power law with 
island area S~i z ,z«0.2 [32, 33, 36]. Other studies 
have found that the density of Internet routers scale as a 
power law with population density, R ~ P a , a « 1.4 [37]. 
In our treatment, the particular exponent is a function 
of the heterogeneity in the Internet architecture deployed 
within a study region. 

In summary, we have presented a theory of firm size 
fluctuations that explains a number of reported statistics 
in economics. The analysis also represents an important 
null hypothesis because it suggests that some emergent 
statistics are a consequence of heterogeneity in the pop- 
ulation of firms. For these systems, understanding and 
quantifying heterogeneity seems to be the central prob- 
lem in understanding macroscopic firm statistics. 
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APPENDIX A: STATISTICS OF BINNING 

Frequently, empirical analysis of data cannot select a 
significant set of firms of a fixed size. The usual remedy 
is to bin by selecting data within a range of sizes. Be- 
cause the range may be large, the conditional statistics 
may change. For example, we examine the conditional 
mean and variance of q, for the subset of firms in the 
bin s m j n < s avg < s max . Within the bin, the distribution 
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of s avg is a complicated function of s n 
denote the average of the log mean size of firms within 
the bin by (s avg \B), the corresponding variance within 
the bin is denoted by var[s avg \B). The conditional ex- 
pectation (q |s avg ) given in Eq. (19) is linear in s avg , con- 
sequently, binning has no impact on the calculation of 
the conditional first moment. The variance of q for the 
set of firms within the bin is given by 

var[<7 \B] = var[q |(s avg \B}} + 4/3 2 var[s avg \B] , (Al) 

where var[g |(s avg \B)] is defined in Eq. (20). The above 
relationship implies that Eq. (23) holds only when we 
have binned firms such that the variance in firms' log 
mean size in each bin, var[s avg \B], is a constant. Fur- 
thermore, because the mean and variance of q is reflected 
in the parameters M and V in Eq. (34), this implies that 
the distribution of growth rates also has a non-trivial de- 
pendence on bin width. 



APPENDIX B: ESTIMATION OF THE GROWTH 
RATE DISTRIBUTION 

The formal expression for the PDF of growth rates for 
a set of firms T follows from Eqs. (31) and (34), 



TUsXF) = / da 



c -[lo g a-Mf/2V r(a+|) 



cosh" 



If) 



In equilibrium, viable firms cannot sustain large sales 
fluctuations, so a > 1 and M 3> V. We replace the 
integrand with its leading-order term in the asymptotic 
expansion in a and perform a change of variables, 



aV2nV 2v^F V2 



/ 



= / dz- 



a z/2-[z-M] 2 /2V 



cosh 



-2 exp(z) 



(Bl) 



We approximate the second integral in Eq. (Bl) using 
the Laplace method to obtain Eq. (35). 

To examine the large-g behavior of Eq. (35), note that 
the Lambert W function is asymptotically logarithmic, 
W(x) ~ log(x). We set w -> \ogg + log V~ + M + V/2 in 
Eq. (35) to obtain, 



n(sl^) 



exp[- log 2 g/2V] 



-jo 



ytjoV + log g 
where the exponent 70 is given by, 



70 = ^ 



1 1 + M + log V 



V 



(B2) 



(B3) 



For finite V, the tail of the distribution of growth rates 
n(<?) decays slower than an exponential, in agreement 
with empirical studies. 
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